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Abstract. £ Oph, HD 149757, is an 09.5 Vnn star with a very high projected rotational velocity (vsini >340 
kms -1 ). It is also a classical runaway star due to its high proper motion. We perform a quantitative analysis of 
its optical spectrum in order to measure important observables of the star such as its mass, effective temperature, 
luminosity and He, C, N, and O abundances. Comparing these observed values to those predicted by the rotating 
evolutionary models of the Geneva group we find that none of the two sets of models is capable of reproducing 
the characteristics of the star. Nevertheless, due to its runaway nature, the reason for this discrepancy may be 
that the star is not the result of the evolution of a single object, but the product of the evolution of a close binary 
system. 
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1. Introduction 

HD 149757, C Oph, is an 09.5 Vnn star (Mason et al lU^5|l . 
with a very high projected rotational velocity (v sin i >340 
kms -1 in all published values, as in Penny. H 9961 Herrero 
et al.,^92 and Howarth & Smith 2001) and a considerable 
He enhancement in its atmosphere, according to Herrero 
et al. (1992) who measured an He abundance by number 
of e(He)=0.16 1 for this object. 

Recent models of massive stellar structure and evolu- 
tion with rotation (Meynet & Maeder, [20001 1217051 Heger 
& Langer, 120001 and references therein) show how stel- 
lar rotation can induce mixing processes in the star that 
bring nuclear processed material from the stellar core to 
the atmosphere, even at the first stages of the life of the 
star, during the main sequence phase (MS). This mix- 
ing is found to be more efficient for more massive objects 
and higher rotational velocities, and it is also a decreas- 
ing function of the metallicity. The sequence of surface 
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contamination in the MS is the following: first a N en- 
hancement appears, that grows as He enrichment starts, 
and finally CO depletions also occur (Meynet & Maeder, 
2000J Heger & Langer. I2000|) . 

Our work (Villamariz et al. 120021 hereafter paper I) 
studies quantitatively the He & CNO abundances for a 
modest sample of 09 stars, showing that surface CNO 
contamination certainly exists for the fast rotator of the 
sample, HD 191423 (vsini=450 kms -1 ). Many other au- 
thors have also found this effect in some OB stars, mainly 
in B supergiants (Gies & Lambert. 119921 Lennon et al., 
[13051 McErlean et al.,HHH3 Daflon et al. l200l"|) . and Venn 
( 1995) also for some A supergiants. But a one-to-one rela- 
tion between stellar rotation and these surface CNO pat- 
terns is still not clear. 

For an object like HD 149757, with such a high rota- 
tional velocity, the He enhancement found in Herrero et 
al. 1)1992(1 agrees also with the above considerations, but 
a study of the CNO abundances is mandatory to give a 
more complete picture. Also the He abundance has to be 
revised to account for the lower He abundances that we 
determine now due to improvements in our analysis (see 
Villamariz & Herrero [5002 and paper I for details). 

Therefore, as the natural progress of our series on the 
chemical composition of Galactic OB stars we present our 
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study about the He & CNO abundances of HD 149757, 
that together with other parameters such as its age and 
mass, will allow us to make a quantitative comparison with 
the predictions of the evolutionary models. 

Nevertheless, the star is a classical runaway (Blaauw 
ITMT1 Gies HWjl due to its high proper motion 

(/i Q =13.07 ± 0.85 mas/yr, ju«j=25.44 ± 0.72 mas/yr, from 
the Hipparcos Catalogue, ESA I199?j> . Therefore, it may 
not be a single star evolving by itself, but the result of the 
evolution of a close binary system, as considered in the lit- 
erature (from Blaauw 1952 to Hoogerwerf et al., 2001). If 
this was the case, our discussion in terms of single massive 
stellar evolution would not be suitable, but the runaway 
origin of ( Oph is not clear in part due to the results of 
this work and also because of another important question. 

The two scenarios accepted nowadays for the forma- 
tion of massive runaway stars, the BSS (Binary-supernova 
scenario) and the DES (Dynamical ejection scenario) as 
explained for example in Hoogerwerf et al. (|2001[1 need 
to be reformulated to account for the recent findings of 
massive stellar evolution with rotation. To explain the ob- 
servable differences between these two types of objects it 
is implicitly considered that surface contamination of the 
runaway can only be due to the accretion of processed ma- 
terial from its primary companion during the evolution of 
the close binary system, which is not necessarily the case 
as explained above: self contamination can occur in the 
isolated evolution of a massive rotating star. 

We will address this question in more detail and at- 
tempt a definitive formation scenario for £ Oph in a forth- 
coming paper; here we restrict our study to the question of 
whether the observable characteristics of C Oph are con- 
sistent with the evolution of a rotating single massive star. 

In sect. ||we present our new analysis of the H/He 
spectrum, that provides the new He abundance together 
with other stellar parameters such as effective tempera- 
ture, logarithmic surface gravity, stellar mass and lumi- 
nosity. In sect. El we determine the CNO abundances and 
present the improvements we have introduced in the anal- 
ysis. In sect. 0] we compare the abundances to what mas- 
sive stellar evolution with rotation predicts for this object 
and finally we conclude in section [3 



2. New analysis of the H/He spectrum 

The improvements that we have introduced in our NLTE 
planeparallcl and hydrostatic analyses of OB stellar atmo- 
spheres, in comparison to those of Herrero et al. ©92), 
are three. First we adopt a slightly different criterion for 
the He I lines suitable for the analysis in the late types 
(09 and early B, see Villamariz 2001 for details). Second 
is the inclusion of line-blocking in the spectral synthesis, 
first as in Herrero et al. I|2000l lb I in table 1) and later on 
with our actual better approximation presented in paper 
I (sect. 3.1, lb II in table 1). The third is the consideration 
of microturbulence (£) in the synthesis of the H/He lines 
(Villamariz & Herrero. HURty . 
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Fig. 1. Evolution of the parameters of HD 149757, from 
those of Herrero et al. ( 1992, model Ml, solid line) to ours 
(model M5, dash-three dotted). Model M2 (dotted) ac- 
counts for the new criterion for Hel lines, M3 (dashed) for 
line-blocking, M4 (dash-dotted) for microturbulence and 
the final M5 (dash-three dotted) accounts for all the previ- 
ous with the updated version of line-blocking. See text for 
details. Our uncertainties are the following: AT o //=1000 
K, A log g=0.10 dex, Ae(He)= 25% and A£=2 kms" 1 



The observations used for this work are those presented 
in Herrero et al. CI992J) , with the ING 2.5m Isaac Newton 
Telescope and the Intermediate Dispersion Spectrograph, 
with the H2400B grating yielding a spectral resolution of 
0.6 A and with exposure times giving a S/N in excess of 
200 in the H^ region. 

In Fig.^we can see the evolution of the stellar param- 
eters of ( Oph from Herrero et al. (1992, model Ml, solid 
line) to the present ones (model M5, dash-three dotted). 
The other three model spectra are the best fits consider- 
ing all the improvements one by one, i.e, first consider- 
ing only the new criterion for the Hel lines (M2, dotted 
line), then adding line-blocking as in Herrero et al. (2000 
M3, dashed), then considering also microturbulence in the 
synthesis of the H/Hc lines (M4, dash-dotted) and finally 
with the updated line-blocking of paper I (M5, dash-three 
dotted). 
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Table 1. Parameters of the models plotted in Fig. Ml to M5, and of all the models used. Note that for CNO e(X) 
is in the standard e(X)=12+log ^, while for He e(He)= jvrm+jv(H e ) • ^ ur vs hii=400 ± 20 kins" 1 . 



Mod. 


T »ff 


logs' 


e(He) 
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e(C) 


e(N) 


e(O) 


Effect 


Line 


Ml 


32 500 


3.55 


0.16 














initial 


solid 


M2 


34 000 


3.70 


0.09 














He I lines 




M3 


34 000 


3.70 


0.09 














lb I 




M4 


34 000 


3.70 


0.09 


15 















M5 


34 000 


3.70 


0.11 


15 











lb II 




M6 


34 000 


3.70 


0.11 


15 


8.19 


7.89 


8.69 






M7 


34 000 


3.70 


0.11 


15 


7.86 


8.34 


8.69 






M8 


34 000 


3.70 


0.09 


15 


7.86 


8.34 


8.69 






M9 


34 000 


3.70 


0.09 


15 


7.86 


8.34 


8.69 







The value of 15 ± 2 kms -1 for microturbulence is the 
mean value found for the two 09 dwarfs analysed in paper 
I, for which microturbulence could be determined indepen- 
dently from their N in and C in lines. 

In all the references quoted at the beginning of the 
section it is shown how the new He abundances deter- 
mined are systematically lower than or equal to the previ- 
ous ones, the greater differences being for the hottest and 
more extended objects. The only exception to this is the 
consideration of the more complete line-blocking of paper 
I, that as can be seen in Table 1 produces a slightly higher 
He abundance. This result is in contradiction to what we 
expected in paper I, where we thought that the new opaci- 
ties would not affect the formation of the H/He lines, since 
the strongest contribution to the H/He occupation num- 
bers comes from a spectral range that is common to both 
treatments of line-blocking. Nevertheles see in Fig. [21 how 
the new treatment produces a (lbll) model with more He I 
and He II ionized into He in. This produces weaker He I and 
He II lines, that explains the necessity to increase the He 
abundance in order to fit the observed lines. 

Regarding the individual and relative quality of the 
line fits, the subsequent models reproduce what we find in 
all previous works. The extremely poor fit to the He n 4686 
line with our final model is perfectly acceptable for this 
line in this spectral range (see discussion on this line in 
Herrero et al.. I2UUU1 sect. 6.2). 

Thus we find that the He abundance of our object is 
nearly normal or slightly enhanced. As clearly shown in 
Table 1 this new abundance relies on the new criterion for 
the He I lines, which is somewhat subjective. Nevertheless, 
this is the only possibility, since the simultaneous fitting of 
all He I lines with a unique set of stellar parameters is not 
possible in this spectral range, as many other authors in 
the literature have also found (see McErlean et al., [1999 
and references therein). 

Howarth & Smith (2001) consider the von Zeipel's vari- 
ation of T e ff and logg with latitude in a rotating object, 
something that is not usually done. Apart from this re- 
markable improvement, their model physics and fitting 
criteria are essentially those of Herrero et al. I| 19921) . and 
therefore their comparable He abundances of 0.20 ± 0.03 
and 0.16 ± 0.04 respectively confirms that our neglection 
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Fig. 2. Comparison of the populations of the ground levels 
of He I (solid), Hen (dotted) and He in (dashed) for the 
two treatments of line-blocking (lbl and lbll, see text). The 
formation depth of all the relevant He lines is above mass 
column 1. 

of the von Zeipel's law in the analysis is not important, 
even for this very fast rotator. 

2.1. Stellar mass, age and position in the HR diagram 

Directly from the previous analysis, making use of the 
obtained parameters and adding to them only the stellar 
absolute visual magnitud (-4.2, Howarth & Prinja 1989), 
we can obtain a stellar radius (as in Kudritzki et al. H980|) 
of 8.3 ± 1.5 R,0, a mass of 19 ± 11 M Q and a stellar 
luminosity of 8.3 ± 4.0 -10 4 L ; we can also place the 
object in the HR diagram to attempt an age determination 
(see Fig. EJ). 

Interestingly, with their more sophisticated method- 
ology, Howarth & Smith l|2001fl find a mass of 20 M Q , 
an area-equivalent radius of 8.5 R© and a luminosity of 
9.12-10 4 L©, in agreement with our results. 

Considering its position in the HR diagram (HRD), 
£ Oph seems to be slightly more evolved than our refer- 
ence slowly rotating 09 V stars in paper I, HD 214680 and 
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Fig. 3. HR diagram with evolutionary tracks of Meynet 
& Maeder lj2U0"Ul dashed and l2UCJ3l solid). The new evolu- 
tionary tracks are for initial equatorial rotational velocities 
of and 300 km s _1 while the older tracks are only for 300 
kms -1 . The position of HD 149757 is shown together with 
objects we studied in paper I (see Table|5|for details), also 
their masses are plotted. The cross at the right shows the 
uncertainties in T c ff and log L. 

Table 2. Parameters of ( Oph and the reference objects 
analysed in paper I 



Star 


Clasif. 


v sin i 


T e {{ 


logs 


e(He) 




HD 149757 


09.5 Vnn 


400 


34000 


3.70 


0.11 


15 


HD 191423 


09 III : n 


450 


35 000 


3.40 


0.12 


20 


HD 214680 


09 V 


50 


37 500 


4.00 


0.10 


15 


HD 34078 


09.5 V 


40 


36 500 


4.05 


0.09 


15 



HD 34078, who lie very close to the zero age MS (ZAMS, 
see Fig. |3J). They are reference objects in the sense that 
they do not show any CNO contamination on their sur- 
face. 

To estimate the age of our object we intersect the ver- 
tical line defined by its effective temperature with the two 
sets of rotating evolutionary tracks (see section [J] for de- 
tails). For the 2000 tracks, there is no intersection with the 
models for 20 M Q and initial rotational velocities of 400, 
500 and 611 kms -1 . These tracks are not shown in Fig. El 
for clarity, but they are the most appropiate considering 
the actual value of vsini > 340 kms -1 of our object. 

Linearly extrapolating the tracks at 200 and 300 
kms -1 for 25 M Q up to 500 kms -1 , the age of our ob- 
ject has an upper limit of 5 Myr (r < 5Myr). Doing the 
same for the 2003 tracks the upper limit shifts to 3.8 Myr. 



3. Analysis of the CNO spectrum 

3.1. Preliminar comparisons 

To address the question of the CNO abundances of our 
target star £ Oph we start by comparing its spectrum to 
that of the fast rotator HD 191423, analysed in paper I. 
Rather than comparing them directly we first degrade the 
spectrum of C Oph, vsini=400 kms -1 , up to the higher 
rotational velocity of HD 191423 (vsini=450 kms -1 , see 
Fig.g. 

In Fig. 0| we see clearly that the C ni-0 n features at 
A 4067-78 of C Oph are stronger than in HD 191423, while 
the Nlll lines at A 4511-24 and A 4379 are weaker. This 
points to less CNO contamination in our target object. 

Since the stellar parameters of both objects are quite 
similar (see Table 2), we can appreciate how different their 
CNO content can be. But of course, as they are not exactly 
the same, the differences in the spectral features will not 
only be due to a different CNO composition. Therefore 
this comparison is just illustrative; only after the full CNO 
analysis can a quantitative comparison be done. 

3.2. Methodology and new He abundance 

The methodology followed to analyse broad-lined spec- 
tra is not the curve of growth method, but direct CNO 
spectral synthesis. Briefly, we compute a full H/He & 
CNO spectrum based on the stellar parameters of our 
object and with the reference CNO abundances of non- 
contaminated objects, to compare to the observed spec- 
trum (model M6). We then compute subsequent model 
spectra until we find the best fit to the observed CNO 
features (for more details see paper I), that in this case 
yields e(C)=7.86±0.3, e(N)=8.34±0.3 and e(O)=8.69±0.3 
in the standard e(X) = 12+log § (model M7, see Fig.EJ). 

The H/He lines are now to be fitted. The previous 
H/He analysis did not consider any metal species for the 
synthesis of the spectra, so the contribution of contiguous 
metallic lines to the H/He lines was neglected. 

For fast rotators however this contributions can be es- 
pecially severe, since the lines are so wide. Therefore, as 
seen in Fig. [HJ the synthesis of the full H/He & CNO spec- 
trum yields a slight correction of the He abundance, from 
0.11 to 0.09. The best H/He model M5 (solid line) turns 
into model M7 (dotted) when the combined H/He & CNO 
spectrum is synthesized (for the CNO abundances pre- 
viously found), and the correction of the He abundance 
to 0.09 produces better fits to the He lines (model M8, 
dashed) . 

The He abundance is determinant in the model atmo- 
sphere, that is held fixed for the subsequent line formation 
of the CNO model spectra. Therefore this change in the 
He abundance, although small, may cause a change in the 
model CNO lines. Therefore we compute a new H/He & 
CNO model spectrum with the previous CNO abundances 
on the basis of the new He abundance, and we find that 
the new fits to the CNO features (model M9, see Fig. 
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Fig. 4. Comparison of HD 149757 (v sini=400 kms" 1 , dotted line) and HD 191423 (v sini=450 kms -1 , solid), with 
the former spectrum degraded to the rotational velocity of the latter. The spectral features used for the CNO analysis 
are indicated, the CNO blend at A 4630-52 is not used as explained in paper I. IB stands for interstellar band. 



are better than before, so the CNO abundances need not 
be changed. These better fits are due also to the improve- 
ments we have implemented in the combined CNO model 
atom used this time for the spectral synthesis (see sect. 

rot . 

3.3. CNO abundances and improvements in our CNO 
model atom 

The CNO abundances of our object are obtained as ex- 
plained in the previous section by fitting the CNO spec- 
tral features found in its spectrum. In Fig. we can 
see the models computed for £ Oph, yielding the follow- 
ing CNO abundances: e(C)=7.86±0.3, e(N)=8.34±0.3 and 
e(O)=8.69±0.3. 

As explained there the determination of the He and 
CNO abundances is coupled. See in Fig. how the small 
change in the He abundance from 0.11 to 0.09 due to the 
consideration of the blends with CNO of the H/He lines, 
does not imply a change in the determined CNO abun- 
dances. 

In this figure we show not only the effect of the He 
abundance on the model CNO lines, but also the effect 
of the new combined CNO model atom used to calculate 
model M9. In paper I we could only determine a lower 
limit for the N abundance of the fast rotator HD 191423, 
because a fit to all the N in features was not possible with 
a unique value of the N abundance. Therefore we have 
improved our combined CNO model atom in the spec- 
tral range around He n 4541 and N ill 4511-24, and around 
He 1 4387 and N III 4379 in order to overcome this difficulty. 

The improvements in each region consisted of includ- 
ing all the Cm, Nm and On transitions of the VALD 



database not already present in the model atom and a re- 
vision of the oscillator strengths for consistency with those 
of VALD (Kupka et al. [13551 Ryabchikova et al. lT555jl . 

In Table 3 we give all the details. The determination 
of the N abundance seems to be resolved, at least for this 
case. 

4. Comparison with evolutionary models 

In the HR diagram in Fig. |3| we see how the evolution- 
ary tracks have evolved from Meynet & Maeder (2000) 
to Meynet & Maeder jg003J . The basic differences be- 
tween the two sets of models are the use of the mass 
loss rates of Vink et al. I|2000l I2001|l and Nugis & Lamers 
(2000) for O and WR stars and the consideration of the 
wind anisotropies induced by rotation, both included in 

Table 3. New transitions included in the CNO model 
atom. The level energies are taken from the NIST database 
(http://physics.nist.gov), in cm -1 , and the loggf values 
from VALD. N in 4379 and Cm 4056 were already in the 
model atom but we have updated their log gf to that of 
VALD 



Trans. 


E ; 


E u 


log gf 


Cm 4056 


324212.49 


348859.99 


0.199 


N m 4379 


320365.7 


343116.4 


1.010 


N m 4527 


287591.5 


309691.8 


-0.471 


N m 4530 


287591.5 


309656.2 


-1.161 


N in 4539 


311690.3 


333712.0 


-0.453 


N m 4544 


311715.2 


333712.0 


-0.152 


N in 4546 


314217.3 


336206.9 


0.004 


N m 4547 


287706.9 


309691.8 


-1.384 
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Fig. 5. H/He & CNO model spectra to determine 
the abundances of £ Oph. First a model with non- 
contaminated CNO composition is computed (M6, solid 
line), and the final best fit to the CNO spectral features 
yields the stellar abundances (model M7, dotted). Both 
models are for e(He)=0.11, £=15 kms -1 and the initial 
CNO model atom (see text). 



the 2003 models. Some other improvements are also intro- 
duced, see Meynet & Maeder 1201131 for details. 

The consequences of these differences in the MS are 
wider tracks at slightly lower luminosities and longer MS 
lifetimes. Concerning the evolution of the rotational veloc- 
ities and surface composition (see Figs. 00 and |3J) the 2003 
models show smoother changes with stellar age: rotation 
decreases considerably more slowly and the He & CNO 
enrichments/depletions appear later and are fewer. 

We compare first our mass and (projected) rotational 
velocity, see Fig. [S] Our mass of 19 M© is in agreement 
with both the 2000 and 2003 values for the models with 
initial 20 M Q mass, although the star's position in the 
HRD lies closer to the 25 M tracks. This effect is related 
to the well-know "mass discrepancy" (Herrero et al. I1992fl 
that is still not completely solved (Repolust et al. 2004). 
Our value of 400 kms" 1 for the projected rotational ve- 
locity is much too high for the 2000 models while a 2003 




Fig. 6. Fitting the H/He lines of £ Oph when their blends 
with CNO lines are considered. The H/He model spectrum 
M5 is shown (solid line) together with the H/He & CNO 
models M7 (dotted) and M8 (dashed), with e(Hc)=0.11 
and e(He)=0.09 respectively. Both M7 and M8 are com- 
puted for the CNO abundances found for this object: 
e(C)=7.86±0.3, e(N)=8.34±0.3 and e(O)=8.69±0.3. 

model starting at w 500 kms -1 would be able to repro- 
duce our measurement. 

Concerning the surface composition (see Fig. [§)), our 
He and O abundances are in agreement with the evolu- 
tionary predictions, that hardly change in the MS, while 
our N enrichment and C depletion are much higher than 
predicted. The comparison is worse with the 2003 models. 

Thus, under the 2003 tracks, our high rotational veloc- 
ity could be explained, but our N enrichment and C deple- 
tion are much too high. The 2000 tracks on the contrary 
cannot explain our high projected rotational velocity, but 
predict our observed abundances of He, N and O. 

This may indicate that either the mixing mechanisms 
are more efficient than considered in the evolutionary 
models for a given rotational velocity, or that the loss of 
angular momentum is less efficient so that high rotational 
velocities last longer in the MS. This is what Herrero & 
Lennon lj2tM|l find for O stars, and Trundle et al. PUU2|| 
for B supergiants. 
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Fig. 7. Fitting the CNO spectral features of ( Oph 
with models M7 (solid line) and M9 (dotted), both 
for e(C)=7.86, e(N)=8.34 and e(0)=8.69. M7 is for 
e(He)=0.11 and the initial CNO model atom and M9 is 
for e(He)=0.09 and the improved one (see text). 

However, until the runaway nature of C Oph is clarified 
(in a forthcoming paper) the results presented here cannot 
serve as constraints for the evolutionary models, since the 
origin of our object may be a close binary system rather 
than a single massive star. 

5. Conclusions 

We find that the fast rotator £ Oph shows a surface com- 
position with N enriched and C depleted, with normal He 
and O abundances. Our abundances are based on our last 
analysis of this object, that includes improvements in our 
method first presented in paper I. In particular we are now 
capable of giving more reliable N abundances due to the 
improvements introduced in the model atom used for the 
analysis. 

When comparing our surface composition with the pre- 
dictions of the evolutionary models of the Geneva group 
(Meynet k Maeder. l2UUUll2UU3|> we find better (although 
not complete) agreement with the 2000 models, while for 



Fig. 8. Comparison of our measured mass and (projected) 
rotational velocity with the MS evolutionary models of 
Meynet & Maeder (|2TjUUl dashed line) and Meynet & 
Maeder (20021 solid). Models are for the quoted pairs of 
initial mass and equatorial rotational velocity. We plot the 
two possible upper limits for the age of our object: 5 and 
3.8 Myrs according to the 2000 and 2003 tracks respec- 
tively. 

the rotational velocity only the 2003 models would pro- 
duce the high value of 400 kms -1 we observe. 

The fact that £ Oph is a runaway prevents us from 
going any further in discussing the meaning of these dis- 
agreements, that may be just indicating that this object 
does not come from the evolution of a single massive star, 
but from the evolution of a close binary system. Only after 
the runaway nature of £ Oph is clarified can this discussion 
continue. 
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Fig. 9. As in Fig0 but for the He & CNO abundances. 
Note that our reference mass fractions, Y(X)o, are those 
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composition in the evolutionary models is solar. 
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